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Abstract: 

We present an algorithm for systematically reconstructing a solution of the {d + 2)- 
dimensional vacuum Einstein equations from a (d + l)-dimensional fluid, extending the non- 
relativistic hydrodynamic expansion of Bredberg et al. in arXiv:1101.2451 to arbitrary order. 
The fluid satisfies equations of motion which are the incompressible Navier-Stokes equations, 
corrected by specific higher-derivative terms. The uniqueness and regularity of this solution is 
established to all orders and explicit results are given for the bulk metric and the stress tensor 
of the dual fiuid through fifth order in the hydrodynamic expansion. We establish the validity 
of a relativistic hydrodynamic description for the dual fluid, which has the unusual property of 
having a vanishing equilibrium energy density. The gravitational results are used to identify 
transport coefficients of the dual fluid, which also obeys an interesting and exact constraint 
on its stress tensor. We propose novel Lagrangian models which realise key properties of the 
holographic fluid. 
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1 Introduction and summary of results 

Black hole physics suggests that gravitational theories are holographic and indeed over the last 
15 years we have witnessed remarkable progress in uncovering holographic relations, mainly 
between [d + 2)-dimensional gravitational theories with negative cosmological constant and 
{d + l)-dimensional local quantum field theories. The case of asymptotically flat spacetimes 
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remains elusive and indeed the structure of the asymptotic solutions [1-3] and the divergences 
of the on-shell action suggest [4] that if there is a dual theory it is nonlocal, see also [5]. 

In a recent paper [6] (see [7-14] for earlier relevant work), a remarkable relation between 
incompressible non-relativistic fluids in (d + 1) dimensions satisfying the Navier-Stokes equa- 
tion and {d + 2)-dimensional Ricci flat metrics was found. The construction of [6] starts 
from considering the portion of Minkowski spacetime between a flat hypersurface Sc, given 
by the equation — = ^rc, and its future horizon 7^"*", the null surface X = T. Then 
the effects of finite perturbations of the extrinsic curvature of Sc, with the intrinsic metric 
7afe kept fixed, were studied and it was found that a regular Ricci flat metric exists provided 
that the Brown- York stress tensor on Sc is that of an incompressible Navier-Stokes fluid. 
More precisely, they work in a hydrodynamic non-relativistic limit and construct the bulk 
metric up to third order in the hydrodynamic expansion. Apart from its intrinsic value, this 
relationship provides a potential example of a holographic duality involving a flat spacetime 
and as such it should be further explored. 

The aim of this paper is two- fold. Firstly, we would like to provide a systematic con- 
struction of the bulk solution to all orders in the hydrodynamic expansion and secondly we 
would like to extract the physical properties of the dual fluid in order to obtain clues about 
the nature of the dual theory. 

Let us summarise our results. First recall that in the hydrodynamic limit essentially any 
{d + l)-dimensional fluid is governed by the incompressible Navier-Stokes equation 

E, = drVi + v^djVi-'nd'^Vi + diP = {), E = dy = o, {i = i,...,d), (i.i) 

where is the velocity field, P is a pressure fluctuation about the background value and ry is 
the kinematic viscosity. This equation has a scaling symmetry and admits a one-parameter 
family of solutions 

vI'\t,x) = evi{e^T,ex), P^'\t,x) = e^P{e'^T,ex). (1.2) 

If one considers now higher-derivative corrections to Ei and E, these are naturally organised 
according to their scaling with e. 

Our first main result is an algorithmic construction of a regular Ricci-flat metric corre- 
sponding to a fluid satisfying (1.1) with specific higher-derivative corrections. We begin with 
the observation that the metric up to order with constant velocity and pressure fields is ac- 
tually flat space in disguise: it can be obtained from the Rindler metric by a linear coordinate 
transformation combining a boost which introduces a constant velocity field with a simple 
shift and rescaling of the radial and time coordinates that introduces a constant pressure 
function P. In fact, these transformations are the most general infinitesimal diffeomorphisms 
that preserve the induced metric at r = Tc and lead to a Brown- York stress tensor which is 
of the form of a fluid in equilibrium. (Recall that the Brown- York stress tensor for Rindler 
spacetime is of that form.) 

To extend the solution to next order in the hydrodynamic expansion, we follow [15] and 
promote the velocity field and the pressure function to be spacetime dependent quantities, 
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subject to the requirement that the Einstein equations hold up to order e^. To satisfy this 
requirement one needs to introduce new terms of order e'^ in the metric and we find that (1.1) 
holds. This reproduces the metric given in [6] (up to a choice of gauge for the dual fluid). 

Next we prove inductively that Einstein equations can in fact be satisfied to arbitrarily 
high order in e, by adding appropriate terms to the metric and modifying the Navier-Stokes 
equation and the incompressibility condition by specific higher derivative corrections. In 
particular, the incompressibility condition E is modified at even powers of e while the Navier- 
Stokes equation Ei is modified at odd powers of e. We carried out this procedure explicitly to 
order and obtained the resultant metric through to order e^, given in (6.5)-(6.6)-(6.8)-(6.9) 
as well as the leading order corrections to E and Ei, given in (6.1) and (6.3), respectively. 

The dual fluid has a number of interesting and unconventional features. Firstly, in equi- 
librium it has zero energy density and positive nonzero pressure. In the hydrodynamic regime, 
the stress tensor receives dissipative corrections, the explicit form of which is given to order 
in (6.10), (6.11) and (6.12). As a consequence of the (bulk) Hamiltonian constraint, the 
fluid stress tensor satisfies the quadratic constraint, 

dTabT^'' = T^. (1.3) 

This is an interesting constraint on the dual theory which will be discussed further below. 

Since the bulk solution is constructed using a non-relativistic hydrodynamic expansion in 
e, the fluid stress tensor appears to be non-relativistic. In turns out however that the complete 
answer, including dissipative terms, can be obtained from a relativistic hydrodynamic stress 
tensor upon using the non-relativistic expansion in e. The rather complicated expressions 
obtained by direct computation actually originate from a much simpler relativistic stress 
tensor. The information encoded in the stress tensor can then be expressed in terms of a 
few transport coefficients. These considerations indicate also that the underlying holographic 
dual theory should be relativistic. 

Recall that a relativistic fluid stress tensor takes the general form 

Tab = pUaUb + phab + H^f,, U^U^f, = 0, (1.4) 

where p and p are the energy density and pressure of the fluid in the local rest frame, is 
the relativistic fluid velocity (normalised such that Uau"- = —1), and hab = lab + UaUb is the 
induced metric on surfaces orthogonal to the fluid velocity. The term 11^^ encodes dissipative 
corrections and may be expanded in gradients of the fluid velocity. When the equilibrium 
energy density is nonzero, one may show that the Landau-Lifshitz prescription [16] implies 
that the energy density p does not receive any gradient corrections. This is in stark contrast 
to the case in hand where the equilibrium fluid has zero energy density and the energy density 
becomes nonzero and negative due to gradient terms. In fact the quadratic constraint (1.3) 
uniquely determines p in terms of the dissipative terms. The leading order result for p is 
given in (7.1). 

From the dissipative terms, II^j^, one can read off the transport coefficients. The rela- 
tivistic corrections are organised according to the number of derivatives they contain. Each 
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of these terms can then be expanded in the non-relativistic hydrodynamic expansion. The e 
expansion may increase the number of derivatives but it never decreases them. This means 
that one can unambiguously extract all transport coefficients but one will have to work to 
sufficiently high order in e in order to obtain all transport coefficients. At first order in deriva- 
tives, there is only one possible term^ and from this the kinematic shear viscosity is obtained. 
If we identify the entropy density of the fluid using the entropy of the Rindler horizon, the 
ratio 77/s takes the universal value, 

T^/s = l/(47r). (1.5) 

At second order, there are six transport coefficients (for a fluid in a flat background) and 
our computation to order determines four of them, see (7.5) and (7.12). The products of 
these coefficients with the (background) pressure are pure numbers and one gets the following 
relations 

2cip = C2P = C3P = Cap = —4. (1.6) 

These results encode all the information contained in the dissipative parts of the stress tensor, 
up to this order. 

Let us now move to the non-dissipative part, 

rrinon— dissipative 7, f-\ rr\ 

T^b =phab- (1-7) 

Note that this satisfies by itself the quadratic constraint. One may ask whether there is a 
Lagrangian model that leads to this stress tensor. It turns out that such a model indeed 
exists and is given by^ 

S = J d'^+^xV^V-W)^- (1-8) 

This describes an ideal fluid with Ua = da(t>/ y—idcj)^ and stress tensor given by (1.7) with 
p = y^—{d(j))'^. The equilibrium conflguration corresponds to the solution <j) = t and near- 
equilibrium configurations are obtained by looking at small fluctuations around this solution. 
The background solution breaks the relativistic invariance and the expansion around this 
solution leads to a non-relativistic nonlocal action (since it contains an infinite number of 
derivatives). Thus, this simple scalar theory (and its generalisations, see (8.13)) reproduces 
the non-dissipative part of the stress tensor and incorporates many of the features ones expects 
from the holographic dual theory. 

The plan of the paper is as follows. In section 2, we discuss a general class of fiat 
metrics admitting a Rindler horizon for which the Brown- York stress tensor on the constant- 
radius slice Sc outside the horizon is that of a perfect fiuid with zero energy density. These 

^For relativistic fluids in general there would also be a bulk viscosity term. For fluids with zero equilibrium 
energy density, this term vanishes as a consequence of the equations of motion (see section 7). 
^Interestingly, this action has also been investigated in the context of dark energy, see [17, 18]. 
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equilibrium backgrounds are used in section 3 as seed metrics to construct corresponding near- 
equilibrium configurations. In section 4, we give the algorithm for systematically constructing 
the bulk solution to arbitrary order in the hydrodynamic expansion, while in section 5, we 
explain the gauge-fixing conditions imposed on the fluid. In section 6, we give explicit results 
for the metric and stress tensor up to order e^, and in section 7, these results are interpreted in 
terms of transport coefficients for the dual fluid. Section 8 proposes a simple dual Lagrangian 
which captures key properties of the fluid, and in section 9, we discuss our results. Finally, 
in the three appendices we present some of the technical results used in the main text. 
In appendix A, we discuss the choice of radial gauge; in appendix B, we derive the class 
of infinitesimal diffeomorphisms preserving the equilibrium form of the Brown- York stress 
tensor; and in appendix C, we derive the general form of the hydrodynamic expansion to 
second order in gradients for a relativistic fluid with vanishing equilibrium energy density. 

2 Equilibrium configurations 

Let us first construct a class of Riemann-flat metrics which have the following properties: (i) 
they admit a co-dimension one hypersurface Sc of flat induced metric 

-fabdx^'dx^ = -rcdr^ + dxidx\ (2.1) 

where the speed of light is arbitrary, x^ = (r, 2:°) and = (r, x*) with i = l,...,d 
(unless otherwise noted, Latin indices are raised with 7*^^ and Greek indices with g^'^); (ii) 
the Brown- York stress tensor [19] on Sc, given by^ 

Tab = ^(^7a6 - Kab), (2.2) 

where Kab is the extrinsic curvature of S^, takes the form of a perfect fluid stress tensor, 

Tab = PUaUb + phab, (2.3) 

and finally; (iii) they are stationary with respect to dj- and homogeneous in the x* directions. 
One such metric with these properties is 

ds^ = g/.^dx^'dx" = -rdr^ + 2drdr + dxidx\ (2.4) 

which describes flat space in ingoing Rindler coordinates. (Upon setting r = 21n(X-|-r) and 
4r = X2 - r2, we recover ds^ — — dT^ -\- d^^.) The induced metric ^afe on the surface 
defined by r = Tc then has the form (2.1), while the Brown- York stress tensor on Sc has the 
form (2.3), where u^- = — r^^, = 0, p = and p = Vc The Rindler horizon is at r = 0. 

We may now obtain further metrics satisfying conditions (i), (ii) and (iii) by considering 
the action of diffeomorphisms on (2.4). Since we wish to obtain a connected thermodynamical 

^From here onwards we will set IGttG = 1. We also assume the validity of classical gravity, i.e., higher-order 
curvature and quantum corrections are suppressed. 
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state space, it is sufficient to consider only diffeomorphisms that are continuously connected 
to the identity, i.e., those obtained by exponentiating infinitesimal diff'eomorphisms. In ap- 
pendix B, we show that there are only two infinitesimal diffeomorphisms yielding metrics 
satisfying conditions (i), (ii) and (iii), after fixing the radial gauge as described in appendix 
A. Exponentiating these, we obtain the following two finite diffeomorphisms. 
The first is a constant boost /3j, given by 



fV^T ^-1^cr-l^ix\ x*^x'-7/3*V^r + (7-l)^x^ (2.5) 



where 7 = (1 — /3^) and /3j = Tc ^^"^vi. The second is a constant linear shift of r and 
associated re-scaling of r, 

r^r-Th, T ^ (1 - ThlvcY^I'^T. (2.6) 

Note that this latter transformation shifts the position of the horizon to r = r/^; as we are 
only interested in the case where Sc is outside the horizon we will restrict r^. < Tc- 

Since these transformations are both linear in the coordinates, one can apply them in 
either order on Rindler spacetime. One obtains 

1 2 dr^ / 2 r-rt \ 27 27^;^ : 

ds = 5— V - J— + . , :drdr , , dx dr 



\- lrc\rc- Th) ' rl{l - v"^ /tc) \l - Th/rc 



+ ( ] dxMr + ( 5., - ;:27r^ ( V— ] ) dxMx^' . (2.7) 



While this metric appears complicated, it is in reality simply fiat space written in a compli- 
cated coordinate system. Nevertheless, the Brown- York stress tensor on Sc has the desired 
form (2.3), where the energy density, pressure and relativistic fiuid velocity are given by 

p = 0, p= J: , J-^ {l,Vi). (2.8) 

In particular, the energy density p vanishes. We stress that there is no fiat solution continu- 
ously connected to Rindler spacetime which generates a nonzero equilibrium energy density 
while at the same time preserving the form of the induced metric on Ec, as shown in appendix 

B. In connection with this, note that the inherent ambiguity in the definition of the Brown- 
York stress tensor [19], which gives the freedom to send Tab Tab + Clab for some constant 

C, produces only a trivial constant shift in the background energy density and pressure given 
by p — 7- p — C and p — ?• p + C . This ambiguity has no dynamical consequences, however, 
and so we will suppress it in the remainder of this paper by setting C = 0. Furthermore, 
since the induced metric ^ab is fiat, there are no additional ambiguities involving the intrinsic 
curvature of Sc. 

Finally, let us consider the thermodynamic properties of the metric (2.7). First, the 
Rindler horizon is a Killing horizon and has a constant Unruh temperature given by 

T= — = — ^ , (2.9) 

27r 47rVrc - rj, ^ ' 
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where the surface gravity n is defined as ^'^V = k^^, where ^ = , ^ {dr + Vidi) is 
normahsed so that ^aS,'^ = — 1 

on X]^. If we then introduce an entropy density associated 
with the Rindler horizon identical to the black hole entropy density s = 1/4G, we note the 
thermodynamic relation^ 

sT = p. (2.10) 

This further implies that two neighbouring equilibrium configurations are related by the 
Gibbs-Duhem relation 

s6T = 5p. (2.11) 

As a final point of interest, we note that the region < r < rh + {l — rh/rc)v'^ is an ergoregion 
with respect to the r-translations, which may lead to effects analogous to superradiance. 



3 Seed metric for near-equilibrium configurations 

Let us now change gears and discuss how to construct near-equilibrium configurations. We 
would like to obtain a solution of the vacuum Einstein equations parameterised by a velocity 
field Vi and a pressure field p which are local slowly varying functions of the coordinates 
while preserving the form of the induced metric on Sc. Working perturbatively in e, in the 
present section we show how to construct this solution up to order e^. In the following section, 
we then show how to extend this 'seed' solution through to all orders in e. 

Clearly, if we promote Vi and p to depend arbitrarily on the coordinates x'^, the metric 
(2.7) is no longer an exact solution of the vacuum Einstein equations. However, treating 
Vi = vf^{T,x) and p = Tc + rc ^^"^ P^^\t,x) as small fluctuations around the background 

— 1/2 

Vi = 0, p = rc in the hydrodynamic limit 



V, 



(^) 



{t,x) = evi{e^T,ex), P^'\t, x) = P{e\, ex), (3.1) 



upon expanding (2.7) (noting that r/j = 2P(^) + 0(e^)) 

we obtain 



ds^ = — rdr^ + 2drdr + dxjdx* 

2(1 Vidx'dT -Ax'dr 



r 

+ ( 1-- 



2 



2P\ 



(u^ + 2P)dr^ + ^dxMx^ +(— + —) drdr + 0(e3), (3.2) 



where each successive line is of one order higher in the e expansion. (Here, and in the following, 
we drop the superscripts on P^''^ and v^f^ for simplicity.) This metric by construction preserves 
the induced metric ^ah on Sc. Moreover, it solves Einstein's equations to O(e^), provided one 



*The relation (2.10) may also be derived by evaluating Komar integrals both at the horizon and at 
Ec. Since the solution is stationary, homogeneous in a;* and Ricci flat, the Komar integrals do not depend 
on the radius of evaluation and one can relate sT defined at the horizon to p defined on using standard 
manipulations, see e.g., [20, 21]. 
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additionally imposes incompressibility, i.e., diV^ = O(e^). The metric (3.2) reproduces the 
solution found in [6] through to order. 

As shown in [6], the Brown- York stress tensor on for the seed solution (3.2) is 

^-1 /y?2 /^7T- r "I 

Tabdx^dx^ = ^-—^ dxMr + ^ dr^ + r^^^ P6ij + ViVj - 2rAvj dxMx^' + 0{e^) . 



(3.3) 

It is interesting to compare this result with the hydrodynamic stress tensor for a relativistic 
fluid with vanishing equilibrium energy density. Up to first order in fluid gradients, this takes 
the form^ 

Tab = PUaUb + phab - ^'qKiab + 0(9^), (3.4) 

where ICab = Ki^t^ic'^d) is the extrinsic curvature of surfaces orthogonal to the fluid velocity^ 
and plays the role of the fluid shear. Expanding this hydrodynamic stress tensor in e to order 
e^, one finds that one can indeed reproduce (3.3) upon setting 

p = + O(e3), p=J- + Z^ + 0{^), r/ = l. (3.5) 

In the following section, we will now describe the algorithm to extend our seed solution 
(3.2) to arbitrarily high order in e, with the pressure and velocity perturbations obeying the 
incompressible Navier-Stokes equations with higher order corrections. Note that, if desired, 
spacetimes describing fluctuations around a background with nonzero velocity and with a 
Rindler horizon located at arbitrary r = may then be found by applying the finite diffeo- 
morphisms (2.5) and (2.6) to the resulting solution. 



4 Constructing the solution to all orders 

Let us begin by assuming we know the bulk metric through order e"~^. The first nonvanishing 
components of the Ricci tensor appear therefore at order e". We now wish to add a new term 
g^iiJ to the metric at order so as to extend our solution to order ©(e""*"^). In addition to 
adding a new term to the metric, we must also consider the effect of gauge transformations 
^("^^ and field redefinitions i^uj"^ and 5P^"'^ also at order e". 
Noting that 

9, ~eO, a, a, ~e2, (4.1) 

adding a new term g\^J at order e" to the bulk metric produces a change in bulk Ricci tensor 

^As noted previously, and as we will discuss in section 7, for a fluid with vanishing equilibrium energy 
density the relativistic divergence K, = dau"^ vanishes at first order in gradients and so the usual bulk viscosity 
term is absent. 

®We reserve the symbol Kab for the extrinsic curvature of Ec itself. 
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at the same order given by 



2 

r 



54-; = -r54^) = - J (a,(r9(:!)) + 25,5^^) - 9,5?^ + 2dlg'^;^^) , (4.2) 
where we write ^j"^ = 6^^ g^J^ and (J-R^f'' = 5'^^5r[J\ 

(n) 

Our goal then is to find a ^/jj/ which cancels out the part of the Ricci tensor at order 
arising from the pre-existing metric up to order e""^. (In fact, as we will see shortly, 
the boundary conditions are such that this requirement leads, after gauge-fixing, to a unique 
solution for gj^) .) The full Ricci tensor at order e" must therefore satisfy 

<^ = ^4"^+^i:^=0, (4.3) 

" (n) 

where RliJ denotes the part of the Ricci tensor at order e" arising from the metric up to order 
e""^. Inspecting the equations for SR^f^J above, we see that such a cancellation will only be 
possible provided the following integrability conditions are satisfied: 

= 5,(4") - ri?W) - = i?W + ri?W. (4.4) 

To verify that these conditions are indeed satisfied, we evaluate the Bianchi identity at 
order e" yielding^ 

o = 5,(4?-r4;!))-i?(;!), 

= a.(4a)+ri?W) ^ 4"J+r4,^) = /i")(r,f). (4.5) 

The integrability conditions are therefore satisfied provided the arbitrary function /i^^ (r, x) 
vanishes. 

Evaluating the Gauss-Codazzi identity on Sc at order e", we find 

V^r,4") = [2V\K^ab-Kai,)]^'''^ = [-2i2,^iV^](") = - J=(i?W+r,4«)) = _ ^ ^(„)(^^ 

(4.6) 

Thus, imposing the conservation of the Brown- York stress tensor on Sc at order e" enforces 
the vanishing of /^"'^(t, x), ensuring the validity of the integrability conditions (4.4). From 
the perspective of the dual fluid, at order (i.e., the first order beyond the seed solution), 
this constraint reduces to the exact Navier-Stokes equation. At subsequent even orders in e. 



^Note that the Ricci tensor 7i^^ is itself of order e" (as we have aheady solved the vacuum Einstein equations 
at order e"~^), so all the remaining covariant derivatives, metrics and inverse metrics appearing in the Bianchi 
identity are those of the background Rindler metric. 
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this step amounts to adding corrections to the incompressibihty condition, while at higher 
odd orders it amounts to adding corrections to the Navier-Stokes equation. 

It is interesting to note that evaluating the Hamiltonian constraint on Sc yields an exact 
and nontrivial constraint on the Brown- York stress tensor of the dual fluid, namely 

dTabT"'' = (4.7) 

Applying this constraint alone to equilibrium fluid configurations, one finds that either the 
equilibrium energy density vanishes peqm = (as is the case here) or else is given by peqm = 
{—2d/{d — l))p, where p is the equilibrium pressure. It would be interesting to see if this 
latter solution branch could be obtained through a modification of the present scenario. We 
will discuss further the implications of the constraint (4.7) for the dual hydrodynamics in 
section 7. 

We may now proceed to use (4.2) to identify a particular solution gj^^J that satisfies the 
Einstein equation (4.3) at order e". One such particular solution is 

~9l^ = 0, 

5^ = /3;"V,x) + (1 - r/r,)/3j")(r,^) + j^'^ dr"(i?Jf) - ri?^ - 2^")), 

5^? = (r, f ) + (1 - r/vM:' ir, ^) " 2 /" dr' dr" b!^-^ , 

J r J r' 

= +Hr/rM':]{r,x) - 2 Tdr'i T dr"R^\ (4.8) 

Here, we have chosen limits so that the integrals vanish for r = rc, facilitating the evaluation 
of boundary conditions to follow shortly. We will leave the lower limit r* appearing in the 

(n) 

inner integral for g^j arbitrary for now. Note however that the rr Einstein equation, in 
combination with the Bianchi identity (4.5), fixes the trace /S^^] = 5*-'73g"j- to be 

pi':l{T,x) = 2r,(i?i") - uRt%^^^. (4.9) 

Allowing now for gauge transformations ^("^'^(r, r, x) at order e", as well as field redefi- 
nitions 6vI"\t,x) and 6P^^\t,x) also at this order, the above solution generalises to 



Jn) 
t/rr 


= 2a^e^")^, 


„H 

Hrr 


rc 


(n) 
9ri 






= ~gt) - e^")'^ + (1 - r/rc)25p(") 


(n) 


= ^(?-2(l-r/r,)<5.f). 


(n) 

4- 





(4.10) 
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Here, to obtain the terms involving Sv^"'^ and (5P^"^ it is sufficient to look at how Vi and 
P appear linearly in the seed solution (3.2), since Svf^"^ ~ (5P(") ~ e". Note that the 
above is still a solution of Einstein's equations, since the additional terms coming from field 
redefinitions and gauge transformations do not contribute to 5R^au , as we see from (4.2). 

(n) 

To fix the gauge, we will impose gru! = 0, for all n > 2. This choice is convenient since 

(n) 

we do not have boundary conditions at Sc for these metric components, unlike for the ff^^ . 
The gr^ components of the full metric are then simply those of the seed solution, i.e., to all 
orders in e, 

P V' 

9rr = 0, g^^ = l + \ , g^.^ = -. (4.11) 

2rc Tc Tc 

Note this gauge choice is fully consistent with our considerations in appendix A (cf. equation 
(A. 2)). To impose this gauge choice then, we must fix 

= _(i (4.12) 

where the ^^"^^(r, x) parameterise the residual gauge freedom. Substituting this back into 
(4.10), we find 

ft¥ = 0, 

g(-) = /3j") - + (1 - r/r,)(/3(") + 5P(")) + j^^ dr' j'^ dr"(i?;f) - rii^!;!) - 2i?('^)), 
= + (1 - -/rc)(/3i? - 26vt^) - 2 r dr' dr"R^-\ 



r 



at^ = /3g + ln(r/r,)/3g - 2 dr'l f dr"E^ . (4.13) 
Imposing the boundary condition = on Sc (so that the induced metric remains 



/3i")(r,x) = ~&>{t,xI /3if{T,x) = 0, /3g(r,x) = 0. (4.14) 



fixed) then fixes 

(n) 

To fix the trace-free part of /3gjj- , we require that the metric must be regular at the future 
horizon r = 0, or equivalently analytic, at each order e". To achieve this, it is useful to choose 
the lower limit in the inner integral for g^J^ in (4.13) to be zero. Since will turn out 
to be a polynomial in r, by setting r^, = we ensure that the outer integral does not generate 
any additional logarithmic terms beyond those associated with the /3g^- ■ (Moreover, the trace 

part of /3g^- is set to zero through (4.9).) Regularity on the future horizon at order e" then 
requires^ that these remaining logarithmic terms vanish: 



/3i-]{r,x) = 0. (4.15) 



A possible caveat here is that we impose order by order regularity but the logarithmic terms could yield 
a finite value on the future horizon when resummed to all orders in e. We thank Bob Wald for discussion on 
this point. 
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In conclusion then, with these boundary conditions, the new part of the bulk metric at 
order e"" is therefore 

gi"^} = 0, 



ffW = (1 - r/r,)F(")(r, x) + T dr' T dr"{R^^ - rR^!:) - 2i?W), 



Tc rr'c 



9- 



(^) = (1 - r/r,)i^(") {t,x)-2 1 dr' / ' dr"Rif 



g^ = -2 TdA r dr"R^, (4.16) 





where the arbitrary functions are 

F(")(T,f) = /3j")(r,x) + 5P(")(r,x), F^''\t,x) = /3i?(T,f) - 25.;^ (r, x) . (4.17) 

We will discuss gauge choices to fix these arbitrary functions in the following section. 

As emphasized above, at each order e" this integration scheme generates a solution that 

" (n) 

is regular on the future horizon r = 0. If Rl^J is a polynomial in r that is regular at r = 0, 
then g^j]) will itself be regular at r = 0. This in turn implies that at the next order R^^^^ will 

" (3) 

be regular as well. Since for the seed solution is regular, this argument ensures regularity 
of gjjii) at r = for each n. 

An additional interesting feature of the hydrodynamic expansion is fact that any vector 
constructed from Vi, P and their derivatives is necessarily of odd order in e, while any scalar 
or rank-two tensor constructed from these variables must necessarily be of even order in e. 
This means that at orders e", where n is odd, the Ricci tensor components R^rr , R^ and 

and the arbitrary function fJ""* (t, x) all vanish identically, and so the only nonzero metric 
component generated by the integration scheme (4.16) is g^-'^ ■ Correspondingly, at orders 
where n is now even, R^^ and Pj^^^ (r, x) are zero, and so and so the only nonzero metric 
components generated by the integration scheme are g^ and glj'^ ■ 

5 Gauge choices for the fluid 

In this section we discuss the gauge choice for the fluid so as fix the arbitrary functions 
Fi"^(r, x) appearing in the integration scheme (4.16). In preparation, let us begin by evalu- 

(n) (n) 

ating the contribution ST^^ to the Brown- York stress tensor arising from the new term gliu 
added to the metric at order e". 

The variation in the extrinsic curvature of Sc at order e" due to gj^) is 



2 — j'^ao 2 ' -'ao 2 
(n) 

where, since 5^^- is already of order e^, only derivatives with respect to r survive, and the 
normal is that associated with the background Rindler spacetime. Thus, we have 

,^,„,_^f^_ ,^,„,_£^^_ ^«<;) = + 1 ,5.2) 
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The new term due to 5^"^ in the Brown- York stress tensor on Sc at order e" is then 



5ri,") = 2(7a.5i^(") - (5.3) 



which evaluates to 



= ^^^^ + ^ /" dr'(<^.,/?l:) - 4")). (5.4) 



rl^'^ V ' c •''0 

The complete Brown- York stress tensor on Sc at order e"" is then 

where T^"^ represents the contribution at order e" due to the metric up to order e"~^. 
5.1 Fixing f/"V> ^) 

From (4.17), we see that the arbitrary vector function F^"'\t,x) appearing at odd orders in 
e is related to redefinitions of the fluid velocity. To fix this ambiguity, we define the fluid 
velocity as the boost from the lab frame to the local rest frame of the fluid, in which the 
momentum density vanishes (i.e., T^-i = in the local rest frame where the fluid velocity is 
purely timelike). In an arbitrary frame where the fluid velocity = 7(1, Wj), this Landau 
gauge condition [16] then reads 

= /l^^^,,u^ = 51 + n'^Uft. (5.6) 

At odd orders in e, the r component of this equation vanishes identically, since there are no 
scalars at odd order. The remaining vector component reads 

o = t(;) + 7;(;-^)^,- + /,("-i)^„ (5.7) 

where p = TahU°'v!' is the energy density in the local rest frame^. Hence at order e" (for odd 
n), we find F-^\t,x) is fixed in terms of known quantities according to 

= t^-^ + f + tI;-^\^ + (5.8) 

where T^^ and p^""^-' denote respectively the full stress tensor and the energy density 
evaluated at order e"~^, both of which are known quantities. 

If we now move to order e""^^ (where n+1 is even), we find the vector component of (5.6) 
vanishes identically (as there are no vectors at even orders), but the r component is nontrivial 
and reads 

= Tj^i)+r^)i;,-r,p("+i). (5.9) 



^While this quantity vanishes for the background solution, it is nonzero at higher order in fluid gradients. 
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Noting that p = j'^{TrT + 2TriVi + TijViVj) and that 1 — rc7^ = — f^7^, one can show that 
(5.9) is equivalent to 

= v'{T(r'^ + T^r'^v,) + r,{T^:^ + T^;'-'\,)v^. (5.10) 

Then, assuming (5.7) is satisfied at order e", (5.9) reduces to 

= t;2(Tjr'^ + T^r^^v, - r,p("-i)). (5.11) 

Thus, (5.9) is satisfied at order e""*"^ provided that it is satisfied at the preceding even order 
e*^"^. Since the seed solution (3.2) explicitly satisfies (5.9) at order e^, we see that the 
gauge condition (5.6) can indeed be consistently imposed at all orders by choosing F^"'\t,x) 
according to (5.8). 

It is quite nontrivial that we can consistently impose the relativistic gauge choice (5.6) on 
the dual fluid, since the latter is constructed according to a non-relativistic expansion scheme 
in which time and spatial derivatives have different order in e. We will return to this topic in 
section 7. 

5.2 Fixing Fj''^(r,f) 

From (4.17), we see that ^^"^(r, x), the arbitrary scalar function appearing at all even orders 
in e greater than two, is related to redefinitions of the pressure fluctuation P. To remove this 
ambiguity, we propose defining the pressure fluctuation P so that the isotropic part of Tij is 
fixed to be 

= + ^) (5.12) 

to all orders. From (5.4), it is clear that we can always choose Fr"\T,x) at every even 
order greater than two so as to ensure that Tij receives no higher order corrections that are 
proportional to 6ij. 

While other gauge choices to define the pressure fluctuation P are possible (see footnote 
12), the present choice is especially simple from an operational point of view. 

6 Results to 0{e^) 

In this section we now apply the integration scheme developed in the preceding sections to 
compute the bulk metric and the ensuing fluid stress tensor at orders up to and including 
in four and five bulk spacetime dimensions. We have checked explicitly at the end of the 
procedure that our solution satisfies the bulk Einstein equations to O(e^) for four and five 
dimensions. 

Since neither the seed metric, nor the algorithm to construct the solution at each step, 
depends on the dimension; and moreover scalars, vectors and rank-two tensors formed from Vi, 
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P have no special properties in three spatial dimensions (unlike in two spatial dimensions^''), 
we expect that the fohowing results are valid in any bulk dimension greater than five as well. 
We have not, however, explicitly checked that the equations of motion hold to O(e^) for bulk 
dimensions greater than five. 

6.1 Corrections to Navier-Stokes and incompressibility 

As noted in [6], imposing conservation of the Brown- York stress tensor at order yields 
the incompressibility condition, while at order we recover the Navier-Stokes equation. In 
section 4, we saw explicitly how conservation of the Brown- York stress tensor at order e" is 
required in order to construct the bulk metric at that same order. 

Applying the integration scheme described in section 4, at order we obtain the following 
corrections to the incompressibility equation: 



1 1 

diVi = —VidiP - Vid'^Vi + -(TijCFij + 0{e^), (6.1) 



where the fluid shear dy- and vorticity cjjj are given by 

(Tjj = 2(9(j-yj) = diVj + djVi, LOij = 'id^iVj^ = diVj - djVi. (6.2) 

As (6.1) is a scalar equation, there are no corrections at odd orders, and so the next corrections 
appear at order e^. 

As expected, the Navier-Stokes equations are modified at order e^: 

drVi + VjdjVi - Tcd'^Vi + diP = --^d'^Vi + 2rcVkd'^dkVi + rcaikdiau ^(^ikdian 

- ^di{akiaki) - ^di{ukiujik) + Vcaudkau - 2vkdkdiP 

- 2{dkVi)dkP - Pd^Vi - ^v'^d'^Vi - ^{dkaii)vkVi + ^{dkUJii)vkVi 
+ 2idkv^)u;kivi + -iP + v^)diP - -VidrP + O(e^) . (6.3) 

Since this is a vector equation, there are no corrections at even orders, hence the next cor- 
rections will be at order e^. 

In writing the form of these correction terms, we eliminated d^P terms using the rela- 
tionship 

d^P+{d,vj){djVi) = 0{e^), (6.4) 

which follows from taking the divergence of Navier-Stokes equation and using incompressibil- 
ity. 



^'^In two spatial dimensions (i.e., d = 2), Vi has two components and dimension dependent identities for rank 
2 tensors exist at order e*. Notably, Oikakj ~ l/25ijakia'ki and u)ikU)kj = l/25ijU)kiL>Jik, where Oij and ujij are 
as defined in (6.2). 
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6.2 Metric at higher orders 

The initial seed metric up to is given in (3.2). Applying our solution algorithm for the 
bulk geometry, at order the only nonvanishing components of the metric are 

9. 



(3) _ {r r,) p^^2 + 2P)^ + Ad^P - (r + r,)d''vi . (6.5) 



2rr 



' c 



Note that only the last of these terms was found in [6]; the first two terms are required in 
order to impose the fluid gauge condition (5.7) at order e^. 

At e^, the only nonvanishing components of the metric are 

g(A = _ (!:_l£)!(^,,^,,) + ^J:_ls^^8v,d\k + ouiau) - t^F(^\ (6.6) 
or^ ore fc 

where ^ 

F,(4) = + ^Pv^ + — - 2r,Vid\i - "^akiTki - 2drP + 2vkdkP, (6.7) 



is fixed by the gauge condition. Also, 

g\j =(1 - -) -2V^Vj{v'^ + 2P) + -V^idj)P - ididjP - -a^kCTkj + -—^UikUJkj 
' c '- ' c ' c ^ ^' c 

+ '^k{i^j)k Vj) + — - — d - -v^idj)V - —Oijiv + 2P) . (6.^ 

At next order we find that the only nonvanishing components of the metric are 



c 



(r + 2rc)(r + 5rc) 4 r - 3rc 2 ^ f , \ a2 , o 
o Vi H — VkO Gik - + rc)VkO Uik + rcaikOiuu 

(r + rc)(r + 4rc) „ 5r + Tc „ , , - 7rcr - 4r^ 

-. ^ikOio-ki — Oi{akiaki) H — di{ujkiujik) 

Arc o lore 

H ^r-^CFkidkCFii -VkdkdiP - aikdkP H i^ikdkP H -d'^ViP 

2 rc rc rc 

+ —V d Vi^ ViVidk(Jki -. dkCTiiVkVi + — dkUJikVkVi + -aikauvi 

2rc rc Arc ^fc 2 

r + 2rc 1 3r — 3rc r + 3rc 

H 7: (^ik^kivi - -iOikCTkivi UikUJkivi H (TkiakiVi 

2rc 2 Atc ore 

_ J. \2 ^ g g 

+ I'^ diV^P - —VkdkPvi - - — P^VkakiP - —drPvi + —aikVkV^ 



ri Vc Tr, rc 

r 2 1 

+ TTo^ikVkV + —VlGklVlVi 

2r^ 2rc 



(6.9) 



Inspecting these results, we see that the bulk metric takes the form of a polynomial in 
r, and as such we would expect it to have only a finite radius of convergence centered about 
r = rc- Note however that this location is arbitrary: by applying the scaling (r, r, Xi,Vi,P) — )• 
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(A^r, T, Xxi, Xvi, }?P) we may set to any chosen value. Let us also mention that, while we 
verified that the d = 4 solution is Petrov II up to order e^^, as noted in [6], the metric fails 
to be of Petrov II type at order e^^. 

Finally, let us comment on the location of the Rindler horizon for the full solution. 
This can be worked out perturbatively in the non-relativistic expansion, as in the analogous 
computation in [22]. Let r = i?(r, Xj) be this position, where R is constructed from fj, 
P and their derivatives. The terms without any derivatives follow from the equilibrium 
solution, which yields Reqm = I'd^ ~ (1 + P/i'c)~'^\- The remaining gradient terms may 
then be determined by imposing at each order in e that the hypersurface be null: = 
g^^d^i{f — R)dy{r — R). To O(e^), there are no gradient terms that can be constructed, hence 
R = 2P + O(e^). To obtain the non-trivial gradient terms appearing at order, however, 
would require knowing g^"^ to order. 

6.3 Stress tensor at higher orders 

The stress tensor of the seed metric was given in (3.3). At order e^, this stress tensor gets 
corrected by the following expression 



2r, 



-3/2 



rcCTikVk - {v"^ + P)vi 



dx'dr . 



(6.10) 



At order e^, the only nonvanishing components of the stress tensor that receive a correction 



are the scalar Trf^ and tensor parts T-J^' as follows. 



V (v + P) - -^CTijaij - rcdijViVj 



dr' 



+ r, 



-5/2 



ViVj{v'^ + P) + 2rcV(^idj)P - ArldidjP - -^(JikOkj 



ik^kj 



3r: 



+ r'tak(^iOJj)k - 2rlv(^i&'^Vj) + ^d^aij - rcV(^idj)V^ - '-^CFijV'^ dxMx^ . (6.11) 



At order e^, we obtain 

.(5) 



Tll>dx<^dx' 



2r-^/2 

2 



-i^Vkd'^aik + ir^VkdkdiP + r^v^d'^Vi + rlviVidkOu + -^cFikdkV^ 

2 

2,, ,, „, I 'c_ „ „, „, „,2/D I „,2\ ^ „,2; 



-^uJikfTkiVi + r'^uJikUJkiVi + -jCTkiakiVi - ViV^{P + v^) - r^v^diP 



rcViVkdkP + rcaikVkV^ + ^uJikVkV^ + ^VkCTkiViVi 



drdx" . 



(6.12) 



To simplify the form of these expressions, we have made use of the constraint equations 
(6.1) and (6.3) below in such a manner that all r derivatives of Vi do not appear in the 
final form of the expressions, as well as equations deriving from these such as (6.4) in such a 
manner that d^P does not appear in the final form of the expressions. 



-17- 



7 Characterising the dual fluid 



Earlier, we noted how the Brown- York stress tensor of the seed metric could be obtained 
from the e-expansion of the hydrodynamic stress tensor for a relativistic fluid. Moreover, 
when constructing the bulk solution, we found one can consistently choose a gauge such that 
u°'Tabh'^ = to all orders in e. Taken together, these observations suggest that it might in 
fact be possible to recover our full gravitational results for the Brown- York stress tensor from 
the e-expansion of some appropriately chosen relativistic hydrodynamic stress tensor. 

Examining the results of the previous section more closely, we find that the energy density 
in the local rest frame is given by 

p = TabU^'u^ = -^^OijUij + 0(e6). (7.1) 

In particular, we see that the energy density vanishes for equilibrium configurations in which 
the fluid velocity is everywhere constant. In the following section, we discuss the theory 
of relativistic hydrodynamics at second order in fluid gradients, paying special attention to 
the modifications necessitated by the vanishing equilibrium energy density. We will then 
proceed to match our proposed relativistic hydrodynamic stress tensor with the Brown- York 
stress tensor derived from our gravitational calculations above, permitting the identification 
of second-order transport coefficients. 

7.1 Relativistic hydrodynamics for vanishing equilibrium energy density 

Defining the relativistic fluid velocity u"" as in (5.6), so that the momentum density in the 
local rest frame vanishes, the fluid stress tensor takes the general form 

Tab = pUaUb + phab + H^fe, u^U^f, = 0. (7.2) 

Here, hab = Jab + UaUb is the induced metric on surfaces orthogonal to the fluid velocity, and 
p represents the pressure of the fluid in the local rest frame. The term 11^^^ encodes dissipative 
corrections and may be expanded in gradients of the fluid velocity. Since in the present case 
the equilibrium energy density vanishes, we must also expand the energy density p in terms 
of gradients of the fluid velocity. Inserting (7.2) into the constraint (4.7) deriving from the 
bulk Hamiltonian constraint, we find 

= p{{d - i)p + 2dp + 2n^) + dni,n^'^^ - (n^)^, (7.3) 

where XI"*" = /i"*n^^. This relation then fully determines the energy density p in terms of 
p and (remembering that we must select the solution branch corresponding to a zero 
equilibrium energy density). Equation (7.3) therefore plays a role somewhat analogous to 
that of the equation of state for a conventional fiuid. 

To write down the gradient expansions for p and Il^^f^ precisely, it is useful to first consider 
the equations of motion at lowest order in fluid gradients: noting that p has no component 
at zeroth order in gradients, these are 

= u''d'Tab = -pdau'' + O{d^), Q = hid'nc Dua = -D^lnp + 0{d^), (7.4) 
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where we have omitted terms of second and higher order in fluid gradients, and we have 
defined D-^ = h'^db and D = u'^da- We may use these equations to simphfy the form of the 
possible coefficients that appear at a given order in the gradient expansion: for example, we 
see that the only possible first order term in the gradient expansion for p, namely dau"", in 
fact vanishes at this order. The expansion for p therefore starts at second order. 

As shown in appendix C, making use of the relations (7.4) we may write down the 
following complete basis for up to second order in gradients, 

nfft = -2r]ICab + ciKlKcb + C2K:ia^\c\b) + c^^a^cb + Cihlh^dcddlnp 

+ cs/Cab D Inp + ceD^ InpD^ Inp + 0{d^), (7.5) 

where t] is the relativistic kinematic viscosity and the ci, C2, etc., are the corresponding 
transport coefficients at second order. Here, we are restricting to flat space^^, and we have 
defined the relativistic shear and vorticity according to 

ICab = Khbd(^^Ud), ^ab = Khbd[cUd], (7-6) 

where symmetrisation and anti-symmetrisation are defined as in (6.2). 

Note also we have not included any terms in (7.5) proportional to some second-order 
scalar times hab itself: this is because our definition of the pressure fluctuation P in (5.12) 
forbids all such terms^^, as may be seen by going to the rest frame in which hij reduces to 

To determine the energy density p in the local rest frame, we now simply insert the 
expansion (7.5) into the constraint (7.3). Expanding to second order in fluid gradients, we 
find 

p = -^/C,5/C'^^ + 0(a3). (7.7) 

p 

7.2 Determination of transport coefficients 

To test the proposed relativistic hydrodynamic expansion given by (7.2), and to identify the 
transport coefficients appearing in (7.5), we now perform a second expansion of (7.5) in e 
up to O(e^). We will then backsubstitute into (7.2), and compare with our results for the 
Brown- York stress tensor in (6.10), (6.11) and (6.12). 

In comparing the two expansions, it is useful to note that the terms associated with the 
coefficients C5 and cg both vanish to O(e^), and so are not constrained by the gravitational 
analysis above. Nevertheless, if required, these transport coefficients may be straightforwardly 
evaluated by extending the gravitational analysis to the appropriate order. 

^^For curved backgrounds additional terms involving the Riemann tensor appear, see e.g., [23]. 

^^Alternative gauge choices for P are however possible. In particular, one may send P ^ P — r^^'^Tl^ /d, 
which would render the dissipative part 11;^^ transverse traceless. This would correspond to the gauge choice 
h'^'T^, = d{r-"'+r-'^'P). 
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In the following analysis, we expect to be able to recover all of the terms in the Brown- 
York stress tensor containing up to two gradients, commensurate with the second-order ac- 
curacy of our corresponding hydrodynamic expansion. Note, however, that the number of 
gradients in a given relativistic term may increase after expanding in e, since we make free use 
of the incompressibility and Navier-Stokes equations to simplify our expressions. Crucially 
though, the order in gradients can never decrease as we go from the relativistic expansion to 
the e-expansion. 

Examining the terms up to order, as noted earlier we find 

P 



1 



P 



+ 



3/2' 



(7.8) 



— 1/2 

i.e., the pressure p of the dual fluid consists of an equilibrium part Vc (as expected from 
the result for Rindler space) and a fluctuation rc ^^^P at order e^. The nature of our gauge 
choice (5.12) is such that this relation is exact and receives no corrections at higher orders. 
The full hydrodynamical stress tensor (7.2) expanded to order is then 



dx 



dr^ + r, 



-3/2 



P5ij + ViVj - 2r]rcdiVj dxMx-^' + O(e^), 

(7.9) 



from which we immediately note rj = 1, although we will leave rj explicit in the following 
formulae. 

Evaluating now the terms at order, we find 



p(3) hydro 
'-ab 



dx'^dx'' = 2r7^/^ 



WcCTikVk - {V^ + P)Vi 



dxMr, 



(7.10) 



which is straightforwardly consistent with (6.10) above. 
At order, we obtain 



rJJ^'^y^^dx'^dx'' 



-3/2 



V {v + P) - rjrcaijViVj - -^(^ij^^ij 



dr' 



+ r, 



-5/2 



HVj{v^ + P) + 27]rcV(^idj)P + C4r^J^didjP + ^r^J^aitakj 



4 



-rjaiji 



4 

dxMx-' . 



(7.11) 



Comparing with Brown- York stress tensor (6.11), we find exact agreement upon setting 

Cl = C2 = C3 = C4 = -4i/r^. (7.12) 

In fact, the only term in (6.11) not captured by the hydrodynamic expansion (7.11) is the 
term proportional to d'^aij: as this term is of third order in fluid gradients, however, we do 
not expect to be able to reproduce it from our second-order hydrodynamic expansion. 
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Finally, at order e^, we obtain 

L o 

-WcViVkdkP + rircaikVkV^ + ^'qu)ikVkv'^ + ^'qvk(TkiViVi drdx". (7.13) 

Checking this result against (6.12), we find that we can indeed reproduce all terms (again, 
apart from the single third-order term proportional to Vkd'^atk) with the assignments (7.12). 

In conclusion then, we have seen how our seemingly complicated results (6.10), (6.11) 
and (6.12) for the Brown- York stress tensor may be recovered from the e-expansion of the 
simple relativistic hydrodynamic stress tensor (7.2), where the relevant transport coefficients 
in (7.5) are given by (7.12). 

8 Models for the dual fluid 

In this section, we present a simple Lagrangian model for the dual fiuid. Since in general one 
would not expect to be able to reproduce the dissipative part of the stress tensor from such 
a model (without first coupling to a heat bath), we focus on reproducing the non-dissipative 
piece of the stress tensor, namely 

Tab=phab. (8.1) 

This corresponds to a fiuid with nonzero pressure but a vanishing energy density in the local 
rest frame, p = Tahu'^vf' = 0. Note that this stress tensor also satisfies the 'equation of 
state' following from the Hamiltonian constraint, dTatT"'^ = T^. The equations of motion 
follow from the conservation of Tab, and are then those in (7.4) but with no higher-derivative 
corrections, i.e., 

d''ua = 0, Dua = -D^\np, (8.2) 

where = /i^Sb and D = u'^da- 

In the previous section, we showed that the stress tensor of the dual fiuid (including 
dissipative terms) may be obtained from a non-relativistic limit of the stress tensor for a 
relativistic fiuid. The fiuid velocity spontaneously breaks relativistic invariance, however, 
and so we are led to consider a relativistic Lagrangian in which the relativistic symmetry is 
spontaneously broken by a background value for the field. The following scalar field action 
satisfies all requirements, 

S = J d'^+^x./^y^^id^. (8.3) 

The field equations are given by 

V'^Ua = 0, Ua = (8.4) 
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where X = —(00)^. Note that Ua satisfies, 



u'^ua = -1. (8.5) 

The stress tensor is given by 

Tab = VXjab + -^da^dbCj) = VXhab, (8.6) 
V A 



This is precisely of the form (8.1) with p = V X and consequently Tab indeed satisfies the zero 
energy density condition, p = TabU°'u^ = 0, and the quadratic constraint dTabT""^ = T^. By 
construction, the stress tensor is conserved and the equations in (8.2) are therefore reproduced. 

Note that (8.4) precisely defines the fluid velocity in terms of a potential (i.e., the fluid 
motion corresponds to potential flow). We also note that taking a generic Lagrangian density 
C{X^ (j)) and then imposing 

= TabU u = (^-2^ + -fabC j ^ = 2X— - C (8.7) 

uniquely picks out the square root action (8.3), up an integration constant which can be 
absorbed into a held redeflnition. 

The equilibrium configuration with pressure p = Vc in the rest frame^^ corresponds to 
the solution, 

= T. (8.8) 

(Recall that the background metric is given by 7afedx"dx^ = — rcdr^ + dxjdx*.) This solution 
spontaneously breaks Lorentz invariance, as required. 

To model a fiuid in its rest frame with small pressure fluctuations about a constant 

— 1/2 

equilibrium value of rc , we must set 

= t + 50(t,x), (8.9) 

The pressure is then given by 

P = ^ = ^{l + 254) + 5^'' - rM,^5(t),if/^ = ^ + (8.10) 

where the last equation serves to define the pressure fluctuation P in terms of the field 
fiuctuation 54), in accordance with (7.8). Similarly, from the components of the relativistic 
fiuid velocity Ur = —rcj = (1 + 54))l\fX and Ui = ^Vi = 54)^il\fX, we may solve for the fluid 
boost velocity Vi in terms of 50, yielding 

- - '^^'^'^ (8.11) 



{1 + 5<P) 



^•^The background solution corresponding to the equilibrium configuration in an arbitrary frame is obtained 
by Lorentz transforming the r.h.s. of (8.8). Note that in such a case Ttt even though the energy density 
is zero. 
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In summary then, we have shown that the non-dissipative part of the fluid stress tensor (8.1), 
describing a fluid with vanishing energy density in the local rest frame, leads naturally to the 
square-root action (8.3). This action is nonlocal in the sense that the expansion around the 
background solution involves an infinite number of derivatives. 

The square root action has various interesting features. The action is real for any timelike 
{d(f>) and the Hamiltonian is positive semi-definite since 

n = X-HdiqyWcP), (8.12) 

so there is no obvious unitarity problem, despite the unconventional nature of the action. 
One can generate a wide variety of field theories by adding extra matter to this scalar field 
model, provided that one expands about zero background values of these fields. For example, 
consider a generic action 

d^+^x^[f[tP, Aa, ^] \/-(a(/>)2 + g[^l;, Aa, 1>] + Aa,n, (8.13) 

with ip denoting fermions, Aa denoting gauge fields and denoting scalars. Here the functions 
f[ip,Aa,^],g[ip,Aa,^] and hlip, Aa,^] are only constrained by the requirement that g vanishes 
with / finite and h either finite or zero when ip = Aa = ^ = 0. Choosing a background in 
which (j) = T with all other fields vanishing will always give a zero energy density fluid, 
as above, for any matter content and choice of functions. Reinstating the l/lGvrG in the 
bulk will introduce a corresponding prefactor in the action above; as usual, a holographic 
correspondence involving classical gravity would be expected to correspond to a dual field 
theory with a large number of degrees of freedom. 

The appearance of the square root in (8.3) might at first sight suggest a connection with 
a brane action. Suppose for example that one considers a ((i-l-l)-dimensional brane embedded 
into a (d -|- 2)-dimensional Minkowski target space. The brane embedding would be described 
by an action 

S = -T [ d'^+^^V^, (8.14) 



where T is the brane tension, ^ are the worldvolume coordinates and = daY'^dbY'^-q^i, 
is the pulled back metric, with the target space coordinates. Fixing static gauge for the 
coordinates (r, x) leads to 

S = -T j d'^+^x^l + (5y)2, (8.15) 

where Y = ig the transverse coordinate to the brane. One can take a tensionless limit 

of this action in which T — )• with = \/TY held fixed to eliminate the constant term in 
the square root 

S = - j d'^+^x^/id^, (8.16) 

but this differs from (8.3) by the absence of the minus term inside the square root. It follows 
that the tensionless limit gives rise to a real action only when {dip)"^ > while the background 
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solution of interest in (8.3) has (dcj))'^ timelike, in which case the onshell brane action would be 
imaginary. One possibility would be to embed the brane into a flat target space of signature 
{d, 2) so that Y ^ iY and ip ^ iip = (p, reproducing (8.3) in the tensionless limit. 

9 Discussion 

In this paper we established a direct relation between (d + 2)-dimensional Ricci-flat metrics 
and (d+l)-dimensional fluids satisfying the incompressible Navier-Stokes equations, corrected 
by specific higher-derivative terms. Our results raise many interesting questions and there 
are diverse directions one may wish to further pursue. 

Perhaps the most interesting question is whether the correspondence extends beyond the 
hydrodynamic regime (on the field theory side) and/or the classical gravitational description 
(on the bulk side) . Is there a string embedding of this correspondence? Even without a string 
embedding, one may ask how this correspondence changes if one adds, for example, a bulk 
stress tensor or considers higher-derivative corrections to Einstein gravity. Will such changes 
modify the properties of the dual fluid? 

The dual fluid has many unconventional features. In particular, it has zero energy density 
in equilibrium but nonzero pressure. What theories have such properties? One of our most 
tantalising observations is that there exists a simple scalar field model that has these properties 
and at the same time has no obvious problems with unitarity, etc. This provides then a 
candidate model for the dual to flat spacetime. Can one obtain this model from branes? Are 
there any other theories which at and near equilibrium are described by such a fluid? 

This paper has focused on the gravity/fluid correspondence, computing transport co- 
efficients of the dual fluid holographically, but it is very exciting to ask how far flat space 
holography can be developed. Defining the holographic theory on Sc, is there a holographic 
dictionary relating bulk computations to quantities computable in the dual field theory model? 

Here, the construction of the bulk metric was achieved using a non-relativistic hydrody- 
namic expansion. We have seen, however, that there is an underlying relativistic description. 
Can one find a manifestly relativistic construction of the bulk metric? 

The metric perturbations are by construction regular at the horizon and throughout 
the region rh < r < Vc, but, in general, one would expect that the series expansion only 
converges in a finite neighbourhood of the surface He- What is the dependence of the radius 
of convergence on the defining data on and what is the range of validity of the gauge choice 
for the radial foliation off the surface Sc? Can one resum the series to obtain a closed-form 
expression for the metric? 

It is important to note here that the hydrodynamic expansion can be made around a 
surface of arbitrary radius Tc, which in particular can be far from the horizon. (As rc — s- oo, 
the pressure and temperature of the equilibrium fluid tends to zero; the stress tensor of the 
fluid is then zero in this limit). The expansion is therefore in general physically distinct from 
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a near-horizon expansion, although it would certainly be useful to understand in more detail 
the limit in which is taken to be close to the horizon, following [7]. 

In this paper, we have used flat space in Rindler coordinates as a seed solution about which 
the hydrodynamic expansion is made. As the present construction uses neither the existence 
of an event horizon, nor the detailed asymptotic structure of flat space, there is no reason 
why it should be applicable only to asymptotically flat or black hole spacetimes. In fact, by 
the equivalence principle, the construction should hold locally in any small neighbourhood. 
Can one patch such a 'local' holographic description of neighbourhoods to obtain a global 
holographic description of general spacetimes? 

Another generalisation would be to start from a different seed solution. This should admit 
at least one Killing vector for the solution to support equilibrium configurations, but the dual 
fluid need not live on a flat spacetime; it could live instead, for example, on a background such 
as Rr X S'^. Note that this would be the geometry of constant-radius slices of Schwarzschild 
spacetime. 

This work also raises many interesting questions from a general relativity viewpoint. We 
have shown explicitly that the data defined on Sc, together with regularity at the horizon, 
suffices to define the solution uniquely as a series expansion in the region between Sc and 
the horizon. Do these solutions exist globally?^^ A priori it would seem far from obvious 
that this radial initial value problem is well-posed and admits a unique solution: it would be 
instructive to understand better these issues. 

All in all, this correspondence provides an interesting arena to explore holography for 
general spacetimes and raises many interesting questions that we hope to return to in the 
future. 

Acknowledgments 

The authors wish to thank Jan de Boer and Michal Heller for discussions. This work is 
part of the research program of the 'Stichting voor Fundamenteel Onderzoek der Materie' 
(FOM), which is financially supported by the 'Nederlandse Organisatie voor Wetenschappelijk 
Onderzoek' (NWO). The authors acknowledge support from NWO; GC and KS via an NWO 
Vici grant, and PM via an NWO Veni grant. 

A Choice of radial gauge 

The boundary conditions imposed on the surface Sc defined by r = Tc do not constrain the 
radial coordinate away from the surface Sc. In the main text, we adopt coordinates such that 
lines of constant form a null geodesic congruence around Sc. Moreover, along each null 
geodesic, we will choose the transverse coordinate r to be an affine parameter. 

^^Note that the same question applies also to the standard AdS/fluid correspondence. In exact parallel, it 
was shown in [15] that the derivative expansion defines the solution uniquely as a series expansion in the region 
between the conformal boundary and the future horizon. It is a non-trivial question whether this provides a 
global solution. We thank Bob Wald for discussions regarding this point. 
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These geometrical conditions can be translated into conditions on the metric components 
as follows. Firstly, the tangent to lines of constant x'^ is n = f{x"',r)dr. Here, n is a null 
vector if and only if g^r = 0. One then has 

n-V.n'^ = nf^drf + /V"a,<7ra. (A.l) 

Hence, vectors n are tangent to geodesies with affine parameter r if and only if drgra = and 
f{x"',r) = /{x""). In summary, we impose the conditions 

g-rr = 0, drgra = 0. (A.2) 

These conditions are similar to those used in the conformal fluid/anti-de Sitter gravity corre- 
spondence (see e.g., [15, 24]). Note that there is in general no guarantee that these coordinates 
may be extended far away from Sc. Rather, these coordinates might break down at some 
finite radius outside Sc if the null congruence has caustics. 



B Infinitesimal diff"eomorphisms preserving equifibrium 

Let us consider a general infinitesimal diffeomorphism around Rindler spacetime. From the 
boundary condition gab = 'Jab at r = rc, we find that on the surface r = rc the following 
conditions hold 

da^j) = 0, (B.l) 
^^e - rcdiC + drC = 0, (B.2) 

-e + 2dr{e-rce) = ^. 

The first of these equations amounts to Killing's equation on d-dimensional Euclidean space 
and is solved by 

= ai{r,T) +h]^ij^{r,T)xK (B.4) 

In particular, we obtain = and djdj^i = 0. Acting on (B.3) with di and using (B.2), 
we obtain 

diC = -2dle. (B.5) 
These equations are integrable for SJ' only if = 0. We then find 

hi] ('^' 't) = 'T'^m (^) + (^) (B-6) 

C = -2dlai{r, t)x' + c(r, r) . (B.7) 
The equation (B.2) is integrable only if djjj] = 0. The equations (B.2)-(B.3) are then solved 

by 
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Finally, the diffeomorphisms preserving the induced metric on Sc are given by 



e = -2d^rai{r,T)x' +c{r,T), 

e = {- -d^Mr,r) + -dra,{r,T))x' + - ^ [ dTc{r,T) + t{r), 

Tc Tc Tc 2rc J 

Ci = ai (r, r) + e^ij] {r)x^ . (B.9) 

One then computes the Brown- York stress tensor on Sc for the metric g^j^ = + C^g^y . The 
resulting stress tensor has the form of a relativistic perfect fluid Tat = pUaUb +p(7af) + UaUb), 
where 

p = 

P= ^ + (-C + 45^0 - 2x'drVi)\r=r, 

Vi = -drtti + 4:dfai\r=r^ ■ (B.IO) 

We have therefore shown that a nonzero energy density cannot be generated from an infinites- 
imal diffeomorphism acting on Rindler spacetime. Note that we did not need to use the gauge 
condition (A. 2) for this argument. 

Let us now require that the fluid is in equilibrium with a uniform velocity profile Vi and 
pressure p on Sc. For convenience, let us define r/j via p = Vc ^^"^{1 + rh/2rc). The set of 
allowed diffeomorphisms then reduces to 

a,(r, r) = Mr) - Vi{r)T + e^'''B,{r) + e-^'''Ci{r), 
c(r, r) = -R{r) + e^/^D{r) + e-^/^E{r), (B.ll) 

where Vi{rc) = Vi and R{rc) = r^- We now fix the gauge freedom associated with redefining 
the surfaces of constant r by imposing the conditions (A. 2). Stationarity and homogeneity 
then further restrict dhgrair, x'^) = 0. The vectors obeying these conditions reduce simply to 
a linear combination of a finite set of transformations, 



2 

^ -^W + 2;^"- — ' 

r = ^io) + ^[ij]^' - + /'^ + e"/'^'' • (B.12) 

Recall that Rindler spacetime can be rewritten in null coordinates u = e"^/^, v = 4re~^/^. 
One then recognises r-translations or alternatively boosts in the {u, v) plane (associated with 
■^(o))' spatial translations (associated with i]^-^), rotations (associated with wjjj]), i)-translations 
(associated with and boosts (associated with i?*), all of which are Killing symmetries 
and do not lead to any perturbation C^g^^. The r-dependent translation (associated with /*) 
amounts solely to shifting gri by If we restrict the gauge choice further so that gri = —Vi/vc 
(as in (4.11)), this then fixes /* = 0. 
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Finally, the nontrivial perturbations are (i) a shift of r combined with a linearised dilata- 
tion of time (associated with r/^); (ii) a linearised boost in the i direction (associated with Vi). 
These two transformations are the linearisation of the two finite transformations (2.5)-(2.6) 

— 1/2 

considered in the main text. The pressure p = rc (1 + fh/'^fc) is the linearisation of the 
pressure (2.8) in the main text. 

C A basis for second-order relativistic hydrodynamics 

In this appendix we show how to construct a general basis for the scalars and transverse 
tensors that appear at second order in the gradient expansion of relativistic hydrodynamics, 
for the special case in which the equilibrium energy density vanishes. We restrict to the case 
in which the dual fluid is in a flat background. 

We begin by writing down all possible symmetric transverse tensors constructed from 
two gradients of the fluid pressure p and velocity u". These take the form X^^ = h'^^^h'^^Xcd 
where the transverse projector h'^ = + u"'Ub, and the possible independent choices for Xcd 
are 

dclnpddlnp, dcddlnp, u'^delnpdcUd, u'^deUcddlnp, 

U^'dedcUd, deU^dcUd, deUcd^Ud, deUcddu", dcUeddu"" (C.l) 

(Note that the term u^dcddUe = —dcU'^ddUe and so is not independent.) We have not con- 
sidered terms consisting of the induced metric hab multiplied by a scalar function of second 
order in gradients since our choice of gauge (5.12) effectively absorbs all such corrections into 
the zeroth order phab term in (7.2). 

We may now reduce this list to just six independent choices by applying the fluid equations 
of motion evaluated to leading order in gradients, given in equation (7.4). In particular, these 
tell us that the fluid divergence /C = dau'^ vanishes at first order in the gradient expansion, 
meaning the term dgU^dcUd in the list above vanishes at second order and may be dropped. 
Using the identity 

daUb = JCab + ^ab " UaDUb, (C.2) 

where the transverse tensors K,ab and Vtab are defined in (7.6) and the fluid derivatives D = 
u'^da and D-^ = h\db, the six remaining independent choices for may be written as 

^a^cb, ^(a^|c|b), ^a^cb, hlhfdcddlap, D^liipD^lup, ICabDlnp. 

(C.3) 

For example, the term h'^^h'l^^u'^ dedcUd reduces to a linear combination of the last three terms 
in the list above, by reversing the order of the derivatives, then pushing the u'^ inside the dc 
and applying (7.4). 

Although in the present case the energy density p in the local rest frame is fixed by the 
constraint (7.3), as a matter of general interest one might wish to construct a complete basis 
for second-order scalar quantities, in a fashion analogous to that done above for symmetric 
transverse tensors. 



- 28 ~ 



One way to construct such a basis is to contract the tensors appearing in (C.l), and 
then apply the leading order equations of motion (7.4) along with the identity (C.2). After 
eliminating all dependencies, we find the following complete set of independent terms: 

/Cab/C"^ ^^afe^^"^ D^lupD^^lnp, DlnpDlnp, D^lnp. (C.4) 

Note that taking the divergence of the leading order equations of motion (7.4) yields the 
result 

rta^ftlnp = D^"" In pD^ Inp - KabK""^ + ^ab^"^ + 0{d^), (C.5) 

which, upon expanding in e, reduces to our earlier result (6.4). 

Using this list of terms, for a general hydrodynamic theory with vanishing equilibrium 
energy density, one could then write 

p = biKabK."^ + b2^ab^''^ + 63-D Inp D Inp + Inp + In p D^" Inp + 0{d^), (C.6) 

where bi, 62, etc., define additional second-order transport coefficients. 
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